t=0:0.1:90;
y=t*pi/180;
U =-(73786976294838206464*atan((3937*(290963886665226059011546923419023514268506438509963372704817829077643468860293120000*sin(4*y) -...
    362246567004494312852946755398053885058044176820794539513106433334067044936581120000*sin(2*y) +...
    2233348235352977769090014541742324563098425698485973534074506825203173031936000000*sin(6*y) +...
    1096312456118655082094780253794332574726338825008092253753644564576833728675840000*sin(8*y)))./...
    (515625*(1885996563306389410393072439710958431408212270682697007743460755587448034178287551360*cos(2*y).^6 +...
    125283908784806749122218958006669528988475227531372979251457362861591733964113149991583*cos(2*y).^4 +...
    297009708733950977754244217387400196605349941608251994122096366475861849200660659024586*cos(2*y).^2 -...
    685439368710424775260689615954375688458115412751135550451020959767670393504090134450409)) - ...
    (3937*(362246567004494312852946755398053885058044176820794539513106433334067044936581120000*sin(2*y) +...
    290963886665226059011546923419023514268506438509963372704817829077643468860293120000*sin(4*y) -...
    2233348235352977769090014541742324563098425698485973534074506825203173031936000000*sin(6*y) + ...
    1096312456118655082094780253794332574726338825008092253753644564576833728675840000*sin(8*y)))./...
    (515625*(1885996563306389410393072439710958431408212270682697007743460755587448034178287551360*cos(2*y).^6 +...
    125283908784806749122218958006669528988475227531372979251457362861591733964113149991583*cos(2*y).^4 +...
    297009708733950977754244217387400196605349941608251994122096366475861849200660659024586*cos(2*y).^2 -...
    685439368710424775260689615954375688458115412751135550451020959767670393504090134450409))).*...
    ((3937*(290963886665226059011546923419023514268506438509963372704817829077643468860293120000*sin(4*y) -...
    362246567004494312852946755398053885058044176820794539513106433334067044936581120000*sin(2*y) +...
    2233348235352977769090014541742324563098425698485973534074506825203173031936000000*sin(6*y) +...
    1096312456118655082094780253794332574726338825008092253753644564576833728675840000*sin(8*y)))./...
    (2578125*(1885996563306389410393072439710958431408212270682697007743460755587448034178287551360*cos(2*y).^6 +...
    125283908784806749122218958006669528988475227531372979251457362861591733964113149991583*cos(2*y).^4 +...
    297009708733950977754244217387400196605349941608251994122096366475861849200660659024586*cos(2*y).^2 - ...
    685439368710424775260689615954375688458115412751135550451020959767670393504090134450409)) - ...
    (3937*(362246567004494312852946755398053885058044176820794539513106433334067044936581120000*sin(2*y) +...
    290963886665226059011546923419023514268506438509963372704817829077643468860293120000*sin(4*y) -...
    2233348235352977769090014541742324563098425698485973534074506825203173031936000000*sin(6*y) +...
    1096312456118655082094780253794332574726338825008092253753644564576833728675840000*sin(8*y)))./...
    (2578125*(1885996563306389410393072439710958431408212270682697007743460755587448034178287551360*cos(2*y).^6 +...
    125283908784806749122218958006669528988475227531372979251457362861591733964113149991583*cos(2*y).^4 +...
    297009708733950977754244217387400196605349941608251994122096366475861849200660659024586*cos(2*y).^2 -...
    685439368710424775260689615954375688458115412751135550451020959767670393504090134450409))).*...
    (4714991408265973525982681099277396078520530676706742519358651888968620085445718878400*cos(4*y).^3 +...
    640564518148831666189042833331179833177937729686985123815362769974864530076902906593115*cos(4*y).^2 +...
    4237081149412375189342579797238529444173813283426369961293613249041441691904075246796890*cos(4*y) -...
    10107555751536686008534272674080886762087912170606619429022810064397862088169184629926005))./...
    (4406882887023435*(843907553957760548443538564097928751471349338147899636572024362832990439865677114245120*cos(4*y).^2 +...
    111357373794524496724705480304009838299109993926604587346790163581747883076455406074068992*cos(4*y) +...
    291430075239833000093061427177536992678984322086769978335436039559514898862865599564873728));
plot(t,U)